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Abstract 

We study the prospects for the measurement of the t - ^ mass difference (AM) and 
the g mass [Mg) in the supersymmetric co- annihilation region at the LHC using tau leptons. 
Recent WMAP measurements of the amount of cold dark matter and previous accelerator 
experiments indicate that the coannihilation region of mSUGRA is characterized by a small 
AM (5-15 GeV). Focusing on taus from X2-^tt— >rrx'j' decays in g and q production, we 
consider inclusive 3r+jet+^ production, with two r's above a high Et threshold and a third 
r above a lower threshold. Two observables, the number of opposite-signed r pairs minus the 
number of like-signed r pairs and the peak of the ditau invariant mass distribution, allow for 
the simultaneous determination of AM and Mg for AM >5 GeV. For example, for AM = 
9 GeV and Mg = 850 GeV and with 30 fb~^ of data, we can measure AM to 15% and Mg to 
6%. 



1 Introduction 



Supersymmetry (SUSY), a symmetry between fermions and bosons, allows for the construction 
of models that link a wide range of physical phenomena. While initially proposed on aesthetic 
grounds, SUSY also allows for cancellation of the quadratic Higgs divergence and, hence, opens 
the window for consistent model building up to the grand unification (GUT) or Planck scales. 
The extension of supersymmetry to a local gauge theory, supergravity [T], led to the develop- 
ment of GUT models [21 El giving a description of physics from the GUT scale down to the 
electroweak scale, in addition to incorporating the successes of the Standard Model (SM). LEP 
data confirmed the vahdity of the idea of SUSY grand unification. Additionally, SUSY models 
with R-parity conservation automatically give rise to a cold dark matter candidate, the lightest 
neutralino, Xv Because these models are consistent to very high energies, they provide descrip- 
tions of the early universe that deeply link particle physics and cosmology. Detailed theoretical 
calculations |1] confirm that GUT models can account for the experimentally observed dark 
matter abundance 

There is good reason to believe that SUSY can be discovered at the Large Hadron Collider 
(LHC) and that it may be possible to determine its parameters with enough precision to 
distinguish among SUSY models that correctly predict the dark matter content of the universe. 
While the ideal is to study the prospects of measuring SUSY in as general a way as possible, 
we consider mSUGRA j2] here, a commonly studied model, for concreteness. In mSUGRA, 
four parameters and one sign determine all the masses and couplings: mo, the universal scalar 
soft breaking mass at Mqut'i "^1/2? the universal gaugino mass at Mqut] ^o? the universal 
cubic soft breaking parameter at Mqut] tan/3= {Hi)/ {H2), where (-^1(2)) is the Higgs vacuum 
expectation value which gives rise to the up (down) quark masses; and the sign of /i, where 
/i is the coefficient of the quadratic term of the superpotential, W = 11H1H2. Dark matter 
abundance and accelerator data constrain the parameter space to only three regions [3j: (1) 
the r- Xi co-annihilation region^ where the Xi is the LSP and AM = M~ — Af^^o is small 

(AM ~ 5-15 GeV), (2) the focus point region where the x\ has a large Higgsino component, 
and (3) the scalar Higgs annihilation funnel region. A bulk region is also allowed by dark matter 
constraints, but it is largely excluded by other constraints. 

A collection of other constraints suggests that, of these regions, the co- annihilation region is 
particularly important. They include the following: the light Higgs and lightest chargino mass 
bounds from LEP [7j, the 6 — > 57 branching ratio bound 8J, and the 3. Sex deviation from the 
SM expectation of the anomalous muon magnetic moment from the muon g — 2 collaboration 
jH]. If the anomalous muon magnetic moment deviation is proven, mostly the co-annihilation 
region will be allowed. Dark matter production and annihilation rates in the early universe are 
believed to control the dark matter relic density we observe today; however, co-annihilation 
can provide an additional mechanism for the reduction of dark matter abundances in the early 
universe. In particular, the r and Xi can co-annihilate sufficiently into SM particles only if 
AM is small. Since, in principle, M~ and can separately range from 100 GeV into the 

TeV domain, if this striking smallness of AM were found at the LHC, it would be a strong 
indirect indication that the x\ is the astronomical dark matter particle. We also note that 

^Note that the notation r here refers to the hghter of the two supersymmetric tau mass eigenstates. 
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many other SUGRA models have a co-annihilation region (since co-annihilation can occur with 
non- universal gaugino masses), so investigating this region has greater applicability than just 
mSUGRA (though, constraints do not necessarily exist in other SUGRA models to select out 
the co-annihilation region). 

The co-annihilation region is difficult to analyze due to the production of low energy r's 
from r —>■ rxi decays. Some previous work has been done to determine how to measure masses 
of SUSY particles in the co-annihilation region. One Monte Carlo study investigated the co- 
annihilation region by fitting various final state mass distributions in order to determine the 
SUSY masses [TUj, but this study did not consider any backgrounds. A second analysis of 
SUGRA models with large tan/5 and large AM (~50 GeV) showed that high Et hadronic r's 
can be used to reconstruct final states at the LHC jTTj. It has also been shown that a small 
AM can be measured at the ILC to 10% ^2], but the LHC will be available long before the 
ILC. An analysis similar to this one has shown that it is possible to measure a small AM at the 
LHC to 18% with a luminosity of 10 fb~^, but it was assumed there that an independent gluino 
mass measurement would exist at the 5% level P^]- While this may be a good assumption, at 
high tan/5 and small AM, low energy final state r's can complicate the interpretation of the 
typical Mg measurment and may affect the measurement uncertainty as well as the reliability 

mi. 

In this letter, we describe new techniques and the prospects of measuring AM and Mg 
simultaneously at the LHC, helping to establish whether we reside in the co-annihilation re- 
gion. Following mi El 5 consider an mSUGRA scenario in the co-annihilation region with 
tan /? = 40 GeV, Aq = 0, and fi > 0, but we allow mo and mi/2 to vary such that AM < 20 GeV. 
In Sec. 2, we describe the requirements to select inclusive 3r-|-jet-|-^ events and discuss the 
differences between this 3r analysis and a similar analysis done with two r's. In particular, we 
motivate our selection criteria and describe the mass and counting observables. In this section, 
it will become clear that a primary experimental requirement for this analysis to be realized 
in practice is the efficiency to identify taus with Et as low as 20 GeV. In Sec. 3, we discuss a 
method to use both observables to simultaneously determine both AM and Mg. Conclusions 
are given in Sec. 4. 

2 Detecting a SUSY Signal in the Co-annihilation Re- 
gion Using the 3t Final State 

The primary SUSY production processes at the LHC are pp—>-qg, qq, gg. In each case they 
decay via q—^q'Xi' or qx2 (or qR—^qXi)'^ 9~^QQ'Xi qqx^] and g—^tti or bbi and their charge 
conjugate states, generally producing high Et jets and gaugino pairs. Since we are interested 
in the t - x\ mass difference, we must investigate the t^tXi decay. The branching ratio of 
Xg— >-rr is about 97% for our parameter space and is dominant even for large mi/2 in the entire 
co-annihilation region. The same is true for the Xi^i^t decay mode. Therefore, the t^tXi 
decay is present in events with XaX2^ xtXi ^ and X2X1 pairs. The decays of xf particles 

are less desirable than the decays of X2 particles even though the decay of squarks and gluinos 
have a larger branching fraction to xf 's than to X2^s. The reason is that X2 decays produce a 
correlated pair of one high energy tau from xH^tt and one low energy tau from t^tXi from 
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which we can extract kinematical information (X2— ^TT^rrx°), whereas decays only produce 
a single low energy tau via x^^tu^tuXi- Therefore, we focus on isolating r pairs from X2 
decays. We note that we only use hadronic r's since the leptonic decays make reconstruction 
and r identification difficult. 

The X2X25 X2Xi^ and X2X1 channels each contain a X2 decay chain and lead to final states 
with four, three, and two r's respectively, each with experimental advantages and disadvantages. 
The 4r final state is the cleanest in that it has very small background, especially from SM 
sources; the disadvantage is that it has poor acceptance and efficiency for reconstructing all 
four r's. The 2r final state is the inverse case; as discussed in there are significant 
backgrounds that must be dealt with, but the acceptance is large because a 2r signature allows 
all three channels. The 3r analysis provides a compromise. We study it here as an alternative 
to the 2r case if the backgrounds for the 2r case are underestimated. It will require actual 
data to determine which analysis performs better or if both are needed. It may be that the two 
analyses can be combined to provide additional information, and this is under study. 

The 3r final state contains both SM and SUSY backgrounds. The primary SM background 
is ti production. This produces two r's, one from each t^Wb-^rub, and a third r coming from 
either a bottom decay (expected to be non-isolated) or from a jet faking a r. There are two 
major sources of SUSY background: real r's from decays of Xi^'s or t^s. Both of these can 
produce low or high energy r's, but they are uncorrelated and, as we will see, can be separated 
statistically from X2 decays by opposite-signed/like-signed subtraction {OS — LS). In addition, 
jets faking r's leads to additional backgrounds due to both SM and SUSY production, but 
these can be handled similarly. 

Our primary method for separating r pairs from X2 decays from other sources is to use 
the standard OS — LS technique, which has been used in other analyses [SI ^M- To do this 
in our case, we order the taus by Et {E^ > E^ > E^) and only consider the pairs tit^ 
and r2r3. For illustration, consider a X2X2 event. We expect ri and r2 to be from the two 
X2^rr decays because the mass difference between the X2 a'^d the r is large (M~o ~267 GeV 

A2 

and M~ ^154 GeV), and we expect to be from the T—>-rXi decay, since the mass difference 

between the r and the Xi is small (M~o ~144 GeV). We would expect a r4 to also come from the 

Xi 

T—^T Xi decay, but we do not consider it in this analysis. Thus, either rirs or r2r3 is from a X2) 
so the taus of one pair are correlated and are opposite-signed (OS). In the other pair the r's are 
uncorrelated and can be OS or LS with equal probability. Therefore, we subtract the number 
of observed LS pairs from the number of OS pairs and only use pairs in the kinematically 
allowed region. As a second example, if our event contains a X2 ^^d a fake r from a jet from 
squark decay, we will have two high energy r's and one low energy r. Again, either the rira 
or r2r3 pair is from the X2, producing a pair that is correlated and OS, while the other pair is 
uncorrelated and has equal probability of being OS or LS. This same method works for any 
background where a X2 is produced along with a r from another source, for example, from xf 
decay. 

To generate our data sets, we simulate our model with all SUSY production using ISAJET 
[in]. We run the generated particles through a detector simulator, PGS [T7], using the CDF 
parameter file for jet finding, and directly use the generated particles for the non-jet objects. 
Finally, a separate Monte Carlo routine simulates the effects of efficiency for r's and fake rate 
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for jets. Based on CDF results and preliminary LHC experimentally studies, we take the r 
efficiency to be 50% for r's with > 20 GeV and the r fake rate for a jet to be 1% ± 0.2% 

Based on the simulations, the full set of selection criteria are listed in Tabled The first two 
r's have high enough average energies that it is reasonable to cut at Et>4:0 GeV and |?7| < 2.5. 
At this energy, we can expect efficient identification. For small AM the energy of the third r, 
from the r^rxi decay, can be very low, often with Et ^ 15 GeV. As in ^S], we loosen the 
Et requirement to 20 GeV as a balance between keeping the maximum number of events while 
still being in a region of reasonable reproduction of the r identification capabilities at the LHC 
detectors. 

To reduce both SUSY and SM backgrounds, we require at least one jet and Since we 
expect all events to end in two x5's, we require the ^ to be large. In addition, as previ- 
ously mentioned, squark decays result in high energy jets. We impose the cuts ^>100 GeV, 
E}^^ ^ > 100 GeV {\r]\ <2.5), and ^+E}^^ ^ > 400 GeV. We require cuts on the jet and ^ 
separately to insure that SM events with unusually high jet Et or fJp have a low likelihood of 
passing cuts. These cuts have the further advantage that events passing them should readily 
pass reasonable jet+^ triggers envisioned for the LHC. 

The invariant mass of the r pair from the X2 decay forms a mass peak and provides excellent 
rejection against both backgrounds fOll^. This can be seen by considering the r pair in the 
chain X2^tt— >rrx^ in the rest frame of the X2- ^^is case, the invariant mass is given by 



where 6 is the opening angle between the r's the kinematic endpoint of this distribution 
corresponds to the 9 = n case. Therefore, the peak and endpoint depends mostly on M~o, 

A2 

M~, and M~o, of which M~o and M~o depend heavily on Mg from the mSUGRA relations. 

Fig. HI shows the M^^^^ distribution for SUSY and ti production. We select r pairs with 
Mrr <100 GeV because this is just beyond the endpoint for all AM and Mg values in our 
parameter range. This makes our results less sensitive to the fake rate uncertainty because jets 
from squarks that fake r's tend to have large Et, which can produce large M^-^-. In Fig. ^ 
we also see that the number of events in the peak of the OS distribution in excess above the 
LS distribution increases with increasing AM. In addition, the peak position also increases 
as a function of AM. Both these features will be used in the next section. We note that the 
ti distribution also has a peak, but the event rate is several orders of magnitude less than the 
SUSY production even without standard isolation cuts, and it is ignored throughout the rest of 
this letter. We also note, that as in Ref. lowering the Et requirement of the to 20 GeV 
is needed for the peak to be visible. 




(1) 
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Invariant tt Mass (GeV) Invariant tt IVIass (GeV) 




Figure 1: The Mtt distribution of OS and LS r pairs for SUSY events with AM = 9 GeV 
(top left) and AM = 20 GeV (top right), and tt background (bottom). The LS distribution 
has a characteristic shape of uncorrelated particles. Note that subtracting the LS distribution 
from the OS distribution leaves a mass peak due to both r's originating from a single x% Foi" 
small AM, this peak contains few counts and is centered at low mass, and for large AM, this 
peak contains many counts centered at a higher mass. There is always a kinematical cutoff 
below ~100 GeV. We note that tt production also produces a peak, but the event rate is several 
orders of magnitude smaller than for SUSY production. Therefore, we assume that the SM 
background is negligible. 
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Table 1: Final selection criteria. 



3 identified r candidates with |?7| < 2.5 and Et > 40, 40 and 20 GeV respectively 
1 jet with Et > 100 GeV and \r]\ < 2.5 
^> 100 GeV 

E^r^' ^ > 400 GeV 

Mrr < 100 GeV where only tiTs and T2T3 invariant mass combinations are considered 



3 Simultaneous Measurement of AM and Mg 

In this section, we define our observables, the number of counts [Nqs-ls) and the ditau mass 
peak position (M^!^^^) and describe their values as a function of both AM and Mg. We then 
show how these two variables can be used to simultaneously measure both AM and Mg, and 
compare our AM result to previous analyses that assume an independent Mg measurement. 

The variable Nqs-ls is the number of LS r pairs subtracted from the number of OS r pairs 
passing all the selection requirements in Table ^ Because we expect the T3 to come from the 
T^TXi decay, the average Et of the r^, its probability of having E^ > 20 GeV, and, therefore, 
Nqs-ls grow with AM. Thus, for a known luminosity, a measurement of Nqs-ls allows for 
a determination of AM. An increase in Mg affects Nqs-ls by decreasing the production rate 
of gluinos, which decreases the number of X2 decay chains produced. However, an increase in 
Mg increases the boost of the which increases the number of accepted X2 decay chains. The 
overall effect is to decrease Nqs-ls with increasing Mg. Though mass changes in squarks and 
the neutralinos also modify both production and boost, these mass effects change with AM 
and Mg in mSUGRA. Therefore, Mg provides a scale for the model. Fig. |2l shows Nqs-ls as a 
function of AM and Mg. We see that Nqs-ls is flat below AM ~5 GeV and nearly linear above 
it. At low AM, the number of r pairs from single decays goes to zero as none of the r's from 
T decay pass the 20 GeV threshold; however, Nqs-ls does not go to zero because of the small 
SUSY background from stop quark pair production and decay via ti— s>(rz/)6x°. 
This background is independent of AM, so the small number of events in the very low AM 
region implies that it is negligible. 

We consider two contributions to the uncertainty in the Nqs-ls measurement: the statistical 
uncertainty and the uncertainty in the fake rate. We note that statistical uncertainty is not the 
usual a/ZV but, from binomial statistics, 2^ jv^°+^ ^" "^^^ ^'-'^ uncertainty in the r fake rate 
also produces a systematic uncertainty on the true value of Nqs-ls- We flnd the uncertainty 
due to the fake rate to be small compared to the statistical uncertainty even though ~20% 
of the Nqs-ls event rate is due to events with at least one jet faking a r. While this may 
seem puzzhng, because of the OS — LS procedure, the additional counts due to the fake rate 
predominantly come from situations where a r pair is produced by a X2 decay, and the flrst 
or second r is due to a jet faking a r. In this situation, the fake rate is bringing r pairs from 
X2 decay that would have been included in the 2r analysis into our sample; therefore, these 
additional counts do not reduce our sensitivity. 
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AM (GeV) Gluino Mass (GeV) 

Figure 2: The left plot shows Nqs-ls as a function of AM with a 1% fake rate with the shaded 
band representing the variation due to the 20% systematic uncertainty on the r fake rate. 
Below ~5 GeV, the third r from t^txI is so soft that there is no signal; therefore, counting 
is dominated by SUSY backgrounds, and the number of counts is flat. Above ~5 GeV, the 
number of counts is nearly linear as a function of AM as more and more r's from t^txI pass 
the 20 GeV threshold. The right plot shows the relationship between Nqs-ls and Mg for AM 
= 9 GeV; Nqs-ls decreases strongly with increasing Mg. 



We define M^^ as the position of the peak of the ditau invariant mass distribution after 
the LS distribution is subtracted from the OS distribution. It directly depends on AM (as 
shown in Fig. and indirectly depends on Mg. The dependences of M^^'^^ enter through Eq. 
1. Because of the dependence of the formula on M~ and M~o, M^^'^^ rises as a function of AM 

(shown in Fig. |2l top left); however, the dependence of M^^^^ on Mg comes indirectly because 
of the mSUGRA relations: Mg ~ 2.8mi/2, M~o ~ 0.8mi/2, and M~o ~ 0.4mi/2. Therefore, 

as Ms changes, so do M~o and M~o, which leads to Mp^^"^ rising as a function of M^. This is 

A2 Al 

shown in Fig. El top right. 

Unlike Nqs-ls the M^!f^^ measurement is dominated by the uncertainty on the number 
of events in the sample until very high luminosities. We take the uncertainty on M^!^'^^ to be 
given by the standard formula, aj^peak = RMS/^/N, where the RMS of the mass distribution is 

TT 

shown in Figs. El bottom left and right, and \/iV is the number of counts in the peak (effectively 
Nqs-ls)j which is a function of the luminosity. For £ = 30 fb~^ and Mg = 850 GeV, we show 
the statistical uncertainty on M^^^"^ as a band in Fig. |21 The uncertainty due to the fake rate 
is unmeasurably small compared to the statistical uncertainty, even for large fake rate values, 
up to 5%. This helps confirm our previous assertion that only r pairs from a X2 survive the 
OS — LS procedure. Pairs not from X2 decays would have an M^^'^^ that is characteristically 
shifted. Since this does not happen, even with high fake rates, we conclude that the additional 
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AM (GeV) Gluino Mass (GeV) 



Figure 3: The upper plots show MP^^'^ function of AM (left) and Mg (right). The lower 
plots show the RMS of the Mrr distribution as a function AM (left) and Mg (right). The 
statistical uncertainty is significantly larger than the uncertainty due to fake rate variation 
and is taken here normalized to expectations for 30 fb~^. Note that the mass method is not 
possible for AM values less than 5 GeV due to the inability to form a mass peak even at high 
luminosities. The M^^^"^ increases strongly with both increasing AM and Mg. 



events are from x^X? or events, which contain r pairs from a X2 ^ind an additional r from 

a fake or a x^. This is confirmed by the Monte Carlo. 

Since Nqs-ls and M"^"^^^ both depend on AM and M^, we can invert these dependences 
to measure AM and Mg simultaneously. To do this, we parameterize Nqs-ls a^^d Mp^^'^ as 
functions of AM and Mg. The contours of constant Nqs-ls and constant Mp^^'^ are shown in 
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Fig.l^for AM = 9 GeV and Mg = 850 GeV. The intersection of these central contours provides 
the measurement of AM and Mg, and the auxihary hues, from expected la uncertainties on 
Nqs-ls and M^^^^, respectively, determine the la region for AM and Mg. Note that we require 
the measurement of Nqs-ls to be statistically significant in order to make a measurement. We 
find that for values of AM > 8 GeV, we require less than 10 fb~^ for 3a significance. For AM 
below 6 GeV, no amount of luminosity suffices to reach the 3a level since the third r is not 
energetic enough to be observed, and no M^^^^ can be constructed. 




Figure 4: The contours of constant Nqs-ls and M^^"-^ for AM = 9 GeV, Mg = 850 GeV, 
and £ = 30 fb~^. The middle lines are the central values while the outer lines show the la 
uncertainty on the measurements. The region defined by the outer four lines indicates the la 
region for the AM and Mg measurements. 

Fig. El shows the expected uncertainty on AM as a function of AM and the expected 
percent uncertainty on AM and Mg as functions of luminosity for AM = 9 GeV and Mg = 
850 GeV. We find that for C = 30 ih-\ we can measure AM to ~15% and Mg to ~6%. It is 
important to note again, however, that our determination of Mg is not a direct measurement, 
but a determination of a parameter in our model, in some sense, the SUSY mass scale of the 
model. It will need to be compared to a direct Mg measurement, assuming one is available. If 
the two results were consistent, it would be a consistency check of the gaugino universality of 
the mSUGRA model and that we are indeed in the co-annihilation region. 

If we assume an independently measured Mg, as in previous analyses f21El, our two ob- 
servables can be considered to be two independent measurements of AM that can be compared 
and combined for further testing of the model. We assume here that Mg=850 GeV and has 
been measured to 5%, and we incorporate this uncertainty into our AM measurement as a 
systematic uncertainty. To do this, we determine the variation in our observables from the ex- 
pected values at a fixed AM but with Mg varying by 5%. We then propagate these variations 
using parameterizations of our observables with Mg fixed as in ^Hj. The results are shown 
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Figure 5: The expected uncertainty on AM and Mg for the simultaneous measurement method. 
The plot on the left shows the uncertainty on AM as a function of AM. The plot on the right 
shows the percent uncertainty in AM and Mg as functions of luminosity for AM = 9 GeV 
and Mg = 850 GeV. At high luminosities, the uncertainties on both AM and Mg continue to 
improve and are statistics limited. 

in the top plots of Fig. IHl We see that for both methods, the systematic uncertainty clearly 
dominates. Combining the two results, using acomUned = / -2 ^ -2 ' ^^"^ ^^^^ with C 



= 30 fb^^ AM can be measured to 12% at our 9 GeV mass point. The mass method is always 
slightly worse than the counting method but does help improve the combined uncertainty. The 
full results are shown in the bottom left and bottom right of Fig. IHl 

4 Conclusions 

We have studied the prospects of simultaneously measuring AM and Mg at the LHC in the 
co-annihilation region using the inclusive 3r + jet + ^ final state for an mSUGRA scenario. 
Because this channel is nearly background free after selection cuts, for AM > 5 GeV we are able 
to form two separate observables: the number of signal events (a X2 cross-section measurement 
equivalent) and M^^'^^. Both vary with AM and Mg allowing the simultaneous measurement 
of both. This is particularly important if we are in the co-annihilation region as it is unclear 
if there will be a high quality Mg measurement available since current methods for measuring 
Mg will likely have to be altered to accommodate low Erp r's in the final state. We find that 
with 30 fb~^ we can measure AM to 15% and Mg to 6% for the example point of AM=9 GeV 
and Mg=850 GeV. While our sensitivity to measuring AM at the LHC is not as good as that 
expected at the ILC, it is quite comparable and should be available much sooner. Further, a 
15% measurement of AM would generally be sufficient to determine if the signal is consistent 
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Figure 6: With the assumption that Mg is measured elsewhere to 5%, we use the Nqs-ls 
and values to make two independent measurements of AM, which can be combined 

to produce a more accurate measurement. The top plots show the sources of uncertainty for 

the counting (left) and the mass (right) methods. In both cases, the systematic uncertainty 
dominates the measurement. The bottom plots show combined results as functions of AM (left) 
and luminosity (right). We note that both methods are systematics limited by £ = 10 fb~^. 



with co-annihilation, and therefore, with the Xi being the dark matter particle. Since no gluino 
mass measurement is possible at the ILC unless a very high energy option is available, this 
may be the only correct measurement of Mg. We also note that we have made no attempt 
to optimize these results, indicating that with actual data from the detector, we will likely 
be able to optimize our cuts, leading to a more precise measurement or a lower luminosity 
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needed for the same sensitivity. We have confirmed that efficient r identification down to an 
Et of 20 GeV is crucial for this analysis as in the 2r analysis ^^l- Further, we expect that 
this analysis and the 2r analysis could be used to complement each other in the establishment 
of a co-annihilation signal at the LHC, and perhaps be combined to produce a more accurate 
measurement. Finally, an analysis of this type can by applied to other SUGRA models provided 
they have a co-annihilation region and do not suppress the production of the X2 particles. 
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